Hybrid stars with color superconductivity within a nonlocal chiral quark model 
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The equation of state for quark matter is derived for a nonlocal, chiral quark model within the 
mean field approximation. Special emphasis is on the occurrence of a diquark condensate which 
signals a phase transition to color superconductivity and its effects on the equation of state. We 
present a fit formula for the Bag pressure, which is densitiy dependent in case when the quark 
matter is color superconducting. We calculate the quark star configurations by solving the Tolman- 
Oppenheimer- Volkoff equations and demonstrate the effects of diquark condensation on the stability 
of hybrid stars for different formfactors of the quark interaction. 
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I. INTRODUCTION 

The investigation of color superconductivity in quark 
matter 0,0 which was revived by applying nonperturba- 
tive QCD motivated interactions 0, finds most of its 
justification in the possible importance for the physics 
of compact star interiors [j| and related observable phe- 
nomena like neutron star cooling [fl 0, gamma-ray 
bursts 0, El 0, IT^ | gr avitational wave signals for com- 
pact stars mergers yj| and others 0|. Since calcula- 
tions of quark pairing predict values of the energy gap 
A ~ 100 MeV and corresponding critical temperatures 
for the phase transition to the superconducting state are 
expected to follow the BCS relation T c = 0.57 A for 
spherical wave pairing, quark matter in compact stars 
should always be in one of the superconducting phases. 

The question arises whether conditions in compact 
stars allow the occurence of quark matter and the for- 
mation of stable configurations of hybrid stars. In order 
to give an answer to this question one has to rely on 
models for the equation of state which necessarily intro- 
duce free parameters and therefore some arbitrariness in 
the results [l^. In particular, the question whether color 
superconductivity shall be realized in the 2-Flavor Super- 
conductivity (2SC) or Color Flavor Locking (CFL) phase 
in the compact star interior has been discussed contro- 
versely |0, ^} . We will restrict our further discussion 
to dynamical models of the NJL type with their parame- 
ters adjusted by fitting hadron properties in the vacuum 
before extrapolating to finite temperatures and densities 
within the Matsubara formalism. 



Within those models it has been shown that strange 
quark matter occurs only at densities well above the 
deconfinement transition, for chemical potentials which 
are barely reached in the very center of a compact star 
EM EU Ell EH ■ The interesting and much investi- 
gated CFL phase could thus play only a marginal role 
for the physics of compact stars. The stability of the so 
obtained hybrid star configurations, however, appears to 
depend sensitively on details of the model, including the 
hadronic phase. 

In the present paper we are investigating this depen- 
dence in a systematic way by employing a nonlocal chi- 
ral quark model which allows to vary the formfactor of 
the interaction kernel while describing the same set of 
hadronic vacuum properties. We provide a polynomial 
fit formula of our quark matter equation of state (EoS) 
which proves useful for applications to compact star phe- 
nomenology, as e.g. the cooling 23] and rotational evo- 
lution 24] or the merging 13] of neutron stars. In order 
to compare the results of the present work for hybrid star 
configurations with observational constraints, we pick the 
example of the compact object RX J185635 — 3754, for 
which limits for both the mass and the radius have been 
reported [25l |26| . A further restriction in the mass - ra- 
dius plane of possible stable configurations comes from 
the constraint given by the surface redshift measurement 
of EXO 0748-676 SI. 



II. EQUATION OF STATE OF HYBRID STAR 
MATTER IN p EQUILIBRIUM 
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A. Quark matter with color superconductivity 

We consider the grand canonical thermodynamic po- 
tential for 2SC quark matter within a nonlocal chiral 
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quark model [TJ where in the mean field approximation 
the mass gap (ft f and the diquark gap A appear as order 
parameters and a decomposition into color (c g {r, &, 3}) 
and flavor (/ 6 degrees of freedom can be made. 

A; { M/c }, T) = £ ^ c,/ (^/> A; M/c T) ,(1) 

c,/ 

where T is the temperature and fi t c the chemical poten- 
tial for the quark with flavor / and color c. 

The contribution of quarks with given color c and flavor 
/ to the thermodynamic potential is 



1 f°° 

2/ d QQ 2 {^[^c{Ef(q)+ /J,f c ),T\ 

n Jo 

u[e c (E f (q) -fi fc ),T]} , 



A 2 



24 G 2 



(2) 



where G\ and G2 are coupling constants in the scalar 
meson and diquark channels, respectively. The disper- 
sion relation for unpaired quarks with dynamical mass 
function rrif(q) — m,f + g(q)<j>f is given by 



Ef(q) = yjq 2 +m? f (q). 
In Eq. (|2J we have introduced the notation 

l + exp(-|) 



lo [e c , T] = Tin 



where the first argument is given by 



e c (0 = cx/t+aitf . 



(3) 



(4) 



(5) 



When we choose the green and blue colors to be paired 
and the red ones to remain unpaired, we have 



A c = g(q)A(8 c , b + 6 Ctg ). 



(6) 



For a homogeneous system in equilibrium, the mini- 
mum of the thermodynamic potential Q q with respect to 
the order parameters {4>f} and A corresponds to a neg- 
ative pressure; therefore the constant £l vac = Y] c ^vac 
is chosen such that the pressure of the physical vacuum 
vanishes. 

The nonlocality of the interaction between the quarks 
in both channels qq and qq is implemented via the same 
formfactor functions g(q) in the momentum space. In 
our calculations we use the Gaussian (G), Lorentzian (L) 
and cutoff (NJL) type formfactors defined as 



9g(q) 
9h(q) 
.9njl(<?) 



exp(-q 2 /A 2 G ) , 
[l + (q/A L ) 2 }-\ 

0(1 - <z/A NJL ) . 



(7) 

(8) 
(9) 



The parameter sets (quark mass m, coupling constant 
G\ , interaction range A) for the above formfactor models 
(see Tab. [|J are fixed by the pion mass = 140 MeV, 
pion decay constant / w = 93 MeV and the constituent 
quark mass m = 330 MeV at T = fi = 32]. The 



diquark coupling constant Gi is a free parameter of the 
approach which we vary as G2 — r) G± . 

Following Ref. [2£| we introduce the quark chemical 
potential for the color c, fi qc and the chemical potential 
of the isospin asymmetry, fij , defined as 



fJ-qc = (Muc + Mdc)/2 

A*/ = (aw - Mdc)/2, 



(10) 
(11) 



where the latter is color independent. 

The diquark condensation in the 2SC phase induces a 
color asymmetry which is proportional to the chemical 
potential fig. Therefore we can write 

AV = H q + -y(<S c ,& + 5c,g ~ 2<5 c , r ) , (12) 

where fi q and fig are conjugate to the quark number den- 
sity and the color charge density, respectively. 

As has been shown in |2j| for the 2SC phase the rela- 
tion 4>u = <t> d = 4> holds so that the quark thermodynamic 
potential is [30| 



Vt q ((f), A; fj, q , fii, fi 8 ,T) + Q vac = — + — 



—o / dqq 2 {c 



£r(-Mg + gA*8 - 



to 



e r(fJ.q - ~A*8 - Pl),T 



e r (-H q + -fig + fij),T 



} 



2 f°° 
2 / dqq 2 {ui 



e b (E(q) - fig - -fig) - fii, T 



e b (E(q) + fi q + -fi 8 ) - fu, T 
e b {E(q) -fi q - ^fig)+fU,T 
e b (E(q) + fi q + -fig) + fij, T 



}, 



(13) 



where the factor 2 in the last integral comes from the 
degeneracy of the blue and green colors (e& = e g ). 

The total thermodynamic potential contains besides 
the quark contribution Vl q also that of the leptons fT d 

n(^,A;/i 9 ,/i/,/i 8) W,T) = Sl q ((f), A; fig, fii, fig,T) 

+ J2 rt d (fn,T) .(U) 

Ze{e,P e ,iA,} 

where latter are assumed to be a massless, ideal Fermi 
gas 



W a {fi,T) 



1 



12tt 2 



7 

180 



7T 2 T 4 . (15) 



At the present stage, we do include only contributions 
of the first family of leptons in the thermodynamic po- 
tential. 
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The conditions for the local extremum of £l q corre- 
spond to coupled gap equations for the two order param- 
eters 6 and A 



dfl 



i=0 o ,A=A o 



dA 



. (16) 



l=<*o,A=A n 



The global minimum of fl q represents the state of ther- 
modynamic equilibrium from which all equations of state 
can be obtained by derivation. 



B. Beta equilibrium, charge and color neutrality 

The stellar matter in equilibrium has to obey the con- 
straints of /3-equilibrium (d — ► u + e~ + i> e , u + e~ — ► 
d+u e ), expressed as 



Mrfc = Hue + Me j 



(17) 



color and electric charge neutrality and baryon number 
conservation. 

We use in the following the electric charge density 



3 ^ 3 

C C 

the baryon number density 



II B 



f,c 



and the color number density 

n 8 = \ ^2(n fb + n fg - 2n fr ) 



(18) 



(19) 



(20) 



The number densities rij occuring on the right hand sides 
of the above Eqs. (fTHf> - (|2*U|> are defined as derivatives of 
the thermodynamic potential l|14|) with respect to corre- 
sponding chemical potentials fij 



dfj,j 



<l> ,A ;T,{ H ,j^i} 



(21) 



Here the index j stands for the particle species. 

In order to express the Gibbs free enthalpy density 
G in terms of those chemical potentials which are con- 
jugate to the conserved densities and to implement the 
j3— equilibrium condition 117|) we make the following al- 
gebraic transformations 



G = ^ ^fc n fc + He-Tie 



f,c 



jjj y^(3/i ge - m){n dc + n uc ) - fi e Q 

C 

1-iB n B + A*Q Q + Ms n s , 



(22) 



where we have defined the chemical potential = 
ij-iq — j-ii conjugate to the baryon number density ub 
in the same way as /iq = — fi e is the chemical potential 



conjugate to Q and /is to ns- Then, the electric and color 
charge neutrality conditions read, 



"8 = 0, 



(23) 
(24) 



at given ns- The solution of the gap equations (|16|) can 
be performed under these constraints. 

The solution of the color neutrality condition shows 
that us is about 5-^7 MeV in the region of relevant 
densities (fj, q ~ 300 -j- 500 MeV) . Since [ii is independent 
of fig we consider \x qc ~ [i q (/ig ~ 0) in our following 
calculations. 

To demonstrate how to define a charge neutral state of 
quark matter in 2SC phase we plot in Fig. Q]the electric 
charge density Q as a function of /jj for different fixed 
values of ^b, when system is in the global minimum of 
the thermodynamic potential and rj — 1. As has been 
shown before in Ref. [23] for the NJL model case, the 
pure phases (A > 0: superconducting, A = 0: normal) in 
general are charged. These branches end at critical values 
of /ii where their pressure is equal and the corresponding 
states are degenerate 

P = PA=o(fJ>B,fJ-I,H e ,T) = PA>o(fJ,B,Hl,He,T) . (25) 

In order to fulfill the charge neutrality condition one can 
construct a homogeneous mixed phase of these states us- 
ing the Gibbs conditions |3lj . 
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FIG. 1: Electric charge density for the 2SC (A > 0) and 
normal (A = 0) quark matter phases as a function of hi for 
different fixed values of \ib- Fhe end points of the lines for 
given fj,B denote states with the same pressure and represent 
subphases in the Glendenning construction. 



The volume fraction that is occupied by the subphase 
with diquark condensation is defined by the charges in 
the subphases 



X — Qa>o/ (Qa>o - Qa=o 



(26) 



and is plotted in Fig. [5]for the different formfactor func- 
tions as a function of /is. 
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In the same way, the number densities for the different 
particle species j and the energy density are given by 



< l 3 = Xn jA>0 + (1 - xKa=o 
£ = X£a>o + (1 - x)£a=o • 



(27) 
(28) 
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FIG. 3: Bag pressure for different formfactors of the quark 
interaction in dependence on the baryon chemical potential 
for r\ = and for r\ = 1. For the latter the superconducting 
phase is realized. 



FIG. 2: Volume fraction \ °f the phase with nonvanishing 
diquark condensate obtained by a Glendenning construction 
of a charge-neutral mixed phase. Results are shown for three 
different formfactors introduced in the text. 



The formulas l|25|l - H28(l define a complete set of ther- 
modynamic relations and can be evaluated numerically. 
In the next section we present the results in a form ana- 
log to the Bag model which has been widely used in the 
phenomenology of quark matter. 



III. QUARK STAR EOS AND FIT FORMULAS 

We calculate the quark matter EoS within this non- 
local chiral model [12| and display the results for the 
pressure in a form remniscent of a bag model 



P (s) = PidM - B^^b) 



(29) 



where P^Ms) is the ideal gas pressure of quarks and 
B^{hb) a density dependent bag pressure, see Fig. [3] 
The occurrence of diquark condensation depends on the 
value of 77 = G2/G1 and the superscript s 6 {5*, TV} in- 
dicates whether we consider the matter in the supercon- 
ducting mixed phase (77 = 1) or in the normal phase 
(77 = 0), respectively. 

According to heuristic expectations, the effect of this 
diquark condensation (formation of quark Cooper pairs) 
on the EoS is similar to the occurrence of bound states 
and corresponds to a negative pressure contribution (Fig. 

EJi. 

For phcnomcnological applications of the quark matter 
EoS (|29|l we provide a polynomial fit of the bag pressure 



b (s Hpb) 



Pid(fJ-B), (J-B < Mc 



(s) 

The coefficients a k as well as the critical chemical 

(s) 

potential fic of the chiral phase transition depend on 
the choice of the formfactor (see Tabs. ITT1 llllll . 

The dependence of the diquark gap on the chemical 
potential can be represented in a similar way by the poly- 
nomial fit 



A(Ms) 



_ J EfcoM^B -A*c ') 



(S)sk 







_ (s) A 6V ) 

Ms < Me 



where the coefficients b k for the different formfactors are 
given in Tab. IIVI 

The volume fraction also can be approximated by poly- 
nomials in the following form 

YLoCk (Ms -Mc X) ) fe , Ms>Mc X) 
X(MS) = { ELo c ^b - vi S) ) k , ^ > Ms > Mi S) (32) 



, [i B < M 



(S) 



The coefficients c k and the chemical potentials fi c are 
given in the Table IVl for different formfactors. 



IV. HADRONIC EQUATION OF STATE AND 
PHASE TRANSITION 

At low densities, quarks will be confined in hadrons 
and an appropriate EoS for dense hadron matter has to 
be chosen. For our discussion of quark-hadron hybrid 
star configurations in the next Section, we use the rel- 
ativistic mean field (RMF) model of asymmmetric nu- 
clear matter including a non-linear scalar field potential 
and the p meson (nonlinear Walecka model), see [34|. 
The quark-hadron phase transition is obtained using the 
Maxwell construction, see Refs. j3f||36] for a discussion. 
The resulting EoS is shown in Fig. 0]for the case 7/ = I 
(left panel) when the quark matter phase is supercon- 
ducting and for 77 = (right panel) when it is normal. 
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FIG. 4: EoS for strongly interacting matter at zero temper- 
ature under compact stars constraints for the coupling pa- 
rameter rj — 1 (left panel) and 77 = (right panel). Dashed 
line: relativistic mean-field model for hadronic matter; solid, 
dash-doted and doted lines correspond to quark matter with 
Gaussian, Lorentzian and NJL formfactor functions, respec- 
tively. 



formfactors and both cases of the diquark coupling, r/ = 1 
(left panel) and 77 = (right panel). 

The EoS of hybrid stellar matter for temperature T = 
is relevant also for calculations of compact star cooling, 
since the star structure is insensitive to the temperature 
evolution for T < 1 MeV. 



V. CONFIGURATIONS OF HYBRID STARS 

In this Section we consider the problem of stability of 
cold (T = 0) hybrid stars with color superconducting 
quark matter core. The star configurations are defined 
from the well known Tolman-Oppenheimer-Volkoff equa- 
tions |37| , written for the hydrodynamical equilibrium of 
a spherically distributed matter fluid in General Relativ- 
ity, see also |34j . 



dP(r) _ [e(r) + P(r)] [m(r) + 4irr 3 P{r)] 



dr 



2m(r)] 



(33) 



where the mass enclosed in a sphere with distance r from 
the center of configurations is defined by 
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m(r) = 4n / e(r')r' 2 dr'. 
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FIG. 5: Baryon number density in units of the nuclear satu- 
ration density as a function of baryon chemical potential fiB- 
Left panel: f] = 1. Right panel: r\ = 0. Line styles correspond 
to the previous figure. 



When comparing the three quark model formfactors 
under consideration, the hardest quark matter EoS is ob- 
tained for the Gaussian, and therefore the critical pres- 
sure and corresponding critical energy densities of the 
deconfinement transition are the smallest, see Table [H] 
The same statement holds for the case 77 = 0, when the 
quark matter phases are normal, see right panel of Fig. 

H 

Acording to the Maxwell construction of the decon- 
finement phase transition, there is a jump in the energy 
density, as is showwn in Fig. 0] 

The corresponding jumps in the baryon densities at the 
critical chemical potentials fig are given in Table IVII 
see also Fig. [S] for the behavoir of ub(ij-b) for all three 



(34) 



These equations are solved for a set of central energy 
densities, see Figs. |S] - 151 An approximate criterion for 
the stability of star configurations is that masses should 
be rising functions of the central energy density e(0). 




500 1000 1500 



E(0) [MeV/fm ] 



FIG. 6: Mass - radius and mass - central energy density re- 
lation for compact stars configurations acording to the EoS 
shown in the left panel of Fig. |1] Hybrid stars with Gaussian 
or Lorentzian quark matter models give stable branches. 



Our calculations show that for the Gaussian and 
Lorentzian formfactors one can have stable configurations 
with a quark core, either with (Fig. EJ) or without (Fig. 
EJ color superconductivity whereas for our parameteriza- 
tion of the NJL model (cutoff formfactor) the configura- 
tions with quark cores are not stable. For the Gaussian 
formfacor case the occurence of color superconductivity 
in quark matter shifts the critical mass of the hybrid star 
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FIG. 7: Same as Fig. |S]for the EoS of the right panel of Fig. 
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FIG. 9: Same as Fig. |S]for Lorentzian formfactor. 



from 1.04 M Q to 0.62 M© and the maximal value of the 
hybrid star mass from 1.83 M Q to 1.79 M . For the 
Lorentzian formfactor the branch of stable hybrid stars 
with 2SC supercondcting quark cores lies in the mass 
range between 1.55 M Q and the maximum mass 1.72 M Q . 
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FIG. 8: Radius-mass relation for Gaussian formfactor includ- 
ing the mass dependence of the quark core radius in both 
cases, one with and the other without 2SC phase. 



ical densities of phase transition to hadronic matter for 
different formfactors of quark interaction. 

We find that the charge neutrality condition requires 
that the quark matter phase consists of a mixture of 2SC 
condensate and normal phase. The volume fraction of 
the condensate phase amounts to 65% - 85% depend- 
ing on the formfactor function of the interaction. In the 
present work we did not consider muons in the quark 
matter phase. Their occurence would increase the vol- 
ume fraction of the superconducting phase by about a 
5%, helping to stabilize the 2SC phase. 

We have shown that for our set of formfactors the NJL 
model gives no stable quark core hybrid stars. The oc- 
curence of the superconducting 2SC phase in quark mat- 
ter supports the stability of the quark matter phase. 

Comparison of the quark core neutron star mass-radius 
relation with the mass and radius of the recently observed 
'small' compact object RX J185635 — 3754 and with the 
constraints from the observation of the surface rcdshift 
for the object EXO 0748-676 shows that our model per- 
fectly obeys those constraints. 

These studies can be viewed as a preparatory step be- 
fore more fundamental nonpcrturbative interactions can 
be provided, e.g. by QCD Schwinger-Dyson Equation 
studies IMlMEl. 



In Figs|S|and[?J|we demonstrate that these models fulfill 
the observational constraints from the isolated neutron 
star RX J185635-3754 [HHl and from the observation 
of the surface redshift for EXO 0748-676 H3- 

The Lorentzian model with normal quark matter has 
marginally stable quark cores with radii less than 2 km, 
in the mass range 1.91 -1.92 M , see Fig. 



VI. CONCLUSION 

We have investigated the influence of the diquark con- 
densation on EoS of quark matter and obtained the crit- 
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TABLE I: Parameter sets (A, Gi A 2 , m) of the nonlocal chi- 
ral quark model for different formfactors discussed in the text. 
The last three columns show the critical temperatures at van- 
ishing chemical potential and the critical chemical potentials 
with and without diquark condensate at vanishing tempera- 
ture, respectively. 
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TABLE II: Coefficients for Bag function fit formula for the 
normal phase case, for different formfactors Eq. 130H . 
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TABLE III: Coefficients for Bag function fit formula for the 
superconducting mixed phase, for different formfactors, Eq. 
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TABLE IV: Coefficients for the diquark condensate fit for- 
mula, for different formfactors, Eq. 13111 . 
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TABLE V: Coefficients for the volume fraction Eq. I32H and 
their valid ranges, for different formfactors. 
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TABLE VI: Limiting densities of the coexistence region be- 
tween quark (Q) and hadron (H) matter phases for differ- 
ent formfactors (no = 0.16 fm -3 is the nuclear saturation 
density) in the first two columns. The third column shows 
the critical baryon chemical potential at the phase transition 
= \ see Fig. |K| The subcolumns indicate the cases 
of superconducting (77 = 1) and normal (77 = 0) quark matter. 



